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METHODS AND APPARATUS FOR DETERMINING A DIRECTION OF ARRIVAL IN 
A WIRELESS COMMUNICATION SYSTEM 



5 FIELD OF THE INVENTION 

The present invention generally relates to a wireless communication system, and 
more particularly to methods and apparatus for determining a direction of arrival (DOA) of 
a signal transmitted by a remote unit in a wireless communication system. 

j 1 0 BACKGROUND OF THE INVENTION 

There are numerous applications that can use the location of a subscriber unit in a 
wireless communication system. For example, the location of the subscriber unit can be 
used for emergency services, so that police, fire and/or medical services can be dispatched 
to the location of the subscriber requesting such an emergency service. In addition, location 

1 5 of the subscriber unit can be used in detecting fraudulent use of the communication system, 
police investigations, and the like. Furthermore, the location of a subscriber unit can be 
used to improve performance of the wireless communication system. For example, wireless 
communication systems that utilize beam-forming transceivers can tune the beams based at 
least in part on the direction of arrival (DOA) of a signal from the subscriber unit. 

20 Wireless communication systems currently have the capability for determining the 

location of a subscriber unit. For example, a subscriber unit can be located within a cell 
through identification of the base station antenna used to provide service in the cell. 
However, as a cell can have a one to three mile radius, the location of the subscriber unit is 
relatively imprecise and impractical for most applications. Additional systems have been 

25 proposed and implemented for determining the location of a subscriber unit, such as the use 
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of a Global Positioning System (GPS) unit at the subscriber unit, triangulation onto a 
transmitting subscriber unit, and Amplitude Difference-based Angle Of Arrival Estimation 
(AD-AOA) as described in U.S. Patent No., 5,786,791, issued to Eugene J. Bruckert on July 
29, 1998 and assigned to Motorola, Inc., which is hereby incorporated by reference. 
5 The AD-AOA estimation provides an accurate location of a subscriber unit in a 

wireless communication system, accounts for multi-path scattering of a received signal and 
unequal receiver gains. The AD-AOA estimation is preferably implemented with a phased 
antenna array having antenna elements that are spaced at a distance that is less than one-half 
of a wavelength of the carrier frequency of the received signal. However, the AD-AOA 
ilO estimation and other location identification techniques and apparatus are preferably 
implemented with any number of phased antenna array configurations, including a 
configuration of a phased antenna array having antenna elements that are spaced at a 
distance that is greater than one-half of a wavelength of the carrier frequency of the received 
signal. 

15 In view of the foregoing, is should be appreciated that it would be desirable to 

provide methods and apparatus for determining a DOA of a subscriber unit in a 
communication system. In addition, it should be appreciated that it would be desirable to 
provide methods and apparatus for determining a DOA of a subscriber unit in a 
communication system with any number of phased antenna array configurations, including a 

20 configuration of a phased antenna array having antenna elements that are spaced at a 
distance that is greater than one-half of a wavelength of the carrier frequency of the received 
signal. Furthermore, additional desirable features will become apparent to one skilled in the 
art from the following detailed description of the drawings and appended claims. 

25 

2 
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BRIEF DESCRIPTION OF THE DRAWINGS 
The present invention will hereinafter be described in conjunction with the following 
drawing figures, wherein like numerals denote like elements, and: 

FIG. 1 is a wireless communication system according to a preferred exemplary 
5 embodiment of the present invention; 

FIG. 2 is the antenna array of FIG. 1 according to a preferred exemplary 
embodiment of the present invention; 

FIG. 3 is an enlarge view of the antenna array of FIG. 2 according to a preferred 
• exemplary embodiment of the present invention; 

;10 FIG. 4 is a side view of the antenna array of FIG. 2 according to a preferred 

exemplary embodiment of the present invention; 

FIG. 5 is a graph illustrating the relationship between the differential electrical phase 
and the direction of arrival (DO A); 

FIG. 6 is a graph illustrating the preferred minimum spacing solution for multiple 
1 5 DO A solutions of various inter-element spacing; 

FIG. 7 is the transceiver with DOA apparatus of FIG. 1 according to a preferred 
exemplary embodiment of the present invention; 

FIG. 8 is an enlarged view of the in-phase (I) and quadrature (Q) components of the 
demodulated signals that are generated by the demodulators of FIG. 7 and provided to the 
20 combiner and DOA computer for subsequent processing according to a preferred exemplary 
embodiment of the present invention; 

FIG. 9 is the DOA computer of FIGs. 7 and 8 according to a preferred exemplary 
embodiment of the present invention; 

FIG. 10 is difference calculator of FIG. 9 according to a preferred exemplary 
25 embodiment of the present invention; and 
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FIG. 1 1 is the angle estimator of FIG. 9 according to a preferred exemplary 
embodiment of the present invention. 
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DETAILED DESCRIPTION OF THE DRAWINGS 
The following detailed description of a preferred embodiment is merely exemplary 
in nature and is not intended to limit the invention or the application and uses of the 

invention. 

r- : 

O FIG 1 is a block diagram of a wireless communication system 20 that includes an 

Q 

K 1 0 apparatus for determining the direction of arrival (DO A) of a remote unit (RU) 22 according 
if; to a preferred exemplary embodiment of the present invention. The wireless 

T communication system 20 described in this detailed description of the drawings is directed 

p 

U to a cellular or personal communication systems (PCS). However, the principles discussed 
■P herein can be readily applied to other wireless based, radio, cable television (CATV), 
15 telephony or satellite telecommunication systems as well as other data communications 
systems. Furthermore, while the communication system 20 of the preferred exemplary 
embodiment of the present invention is configured to utilize a Code Division Multiple 
Access (CDMA) protocol as described in TIA/EIA Interim Standard IS-95A, Mobile 
Station-Base Station Compatibility Standards for Dual-Mode Wideband Spread Spectrum 
20 Cellular Systems, Telecommunications Industry Association, Washington, D.C. July 1993 
(hereinafter referred to as the "IS-95A Reference Document"), which is hereby incorporated 
by reference, the communication system 20 may utilize other system protocols such as, but 
not limited to, the Narrowband Advanced Mobile Phone Service (NAMPS) protocol, the 
Advanced Mobile Phone Service (AMPS) protocol, the Personal Digital Cellular (PDC) 
25 protocol, the Global System for Mobile Communications (GSM) protocol, Two-Way Paging 

4 
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protocols, or the United States Digital Cellular (USDC) protocol and future generations 
thereof. 

In addition to the remote unit 22, the communication system 20 includes, but is not 
limited to, a sectorized base site 24 having a phased antenna array 26 and transceiver with a 
5 DOA apparatus 28 in operative communication with another communication system, such 
as a public switched telephone network (PSTN) 30. The remote unit 22 transmits a radio 
frequency (RF) signal 32 that can be received by a first antenna 34, second antenna 36, third 
antenna 38, fourth antenna 40, fifth antenna 42 and sixth antenna 44 of the antenna array 26. 
p The antennas (34,36,38,40,42,44) of the base site 24 can be configured to form the phased 
IjlO antenna array 26, or the antennas (34,36,38,40,42,44) of the base site 24 can be used in 
jjj conjunction with an antenna or multiple antennas of other base sites to form the antenna 
array 26. The RF signal 32 received by the phased array antenna 26 is provided to the 
L ; transceiver with DOA apparatus 28, which decodes and demodulates the RF signal 32 and 
% identifies the DOA of the RF signal 32. 

' W 15 Referring to FIG. 2, FIG. 3 and FIG. 4, the antenna array 26 is configured to have a 

first side 46, a second side 48 and a third side 50 that provide a triangular supporting 
structure for the antennas (34,36,38,40,42,44). The antennas (34,36,38,40,42,44) are 
preferably placed at or near the corners of the three sides (46,48,50), and antenna pairs at or 
near the corners (47,49,51) (e.g., the first antenna 34 and the second antenna 36, the third 
20 antenna 38 and fourth antenna 40, and the fifth antenna 42 and the sixth antenna 44) are 
used to identify a DOA estimate of the signal transmitted by the remote unit. However, any 
number of antennas, antenna array configurations, antenna pairings, sides and supporting 
structure shapes can be utilized in accordance with the present invention. 

There are two sets of antenna array characteristics that are preferably used for DOA 
25 estimation in accordance with the preferred embodiment of the present invention. The first 

5 
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antenna array characteristic is the positioning of the individual antenna elements within the 
array. The positioning of the individual antenna elements within the array yields phase 
differences between the signals received on each of the antennas, which are a function of the 
antenna positions and the DOA. (See, Hamid Krim and Mats Viberg, Two Decades of 
5 Array Signal Processing Research, IEEE Signal Processing Magazine, Vol. 13, No.4, July 
1996, for a description of the relationship between phase differences and DOA, which is 
hereby incorporated by reference.) The vector-valued collection of the phase differences 
can be expressed as: 

|j a(G)=[al(G)a2(0)...aN(9)]T (1) 

ffllO Where 9 is the DOA, ai(9)=exp(j(pi(9)), cpi(B) is the complex valued phase shift of antenna 
W element i, and j=sqrt(-l). 

The second antenna array characteristic that enables the estimation of the DOA is the 

Q 

n\ antenna directivity (i.e., the patterns of the antenna elements of the array). Each antenna 
O element of the array has a gain and phase shift associated with each DOA. The vector- 
15 valued collection of the antenna element directivities for an antenna array can be expressed 
as: 

b(9)=[b 1 (9) b2(9) . . . bN(0)]T (2) 
Where bi(0) is the directivity pattern of antenna element i. The directivity patterns of the 
antenna elements can be modified from the isolated value for each antenna element using a 
20 mutual coupling configuration between the elements. This mutual coupling configuration 
between the elements can be utilized in accordance with the present invention with the use 
of directivities measured for the antenna elements while positioned in the desired array 
geometry. 

The combination of the two characteristics due to the antenna element positions and 
25 directivities forms a vector-valued response function commonly referred to as the 

6 
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composite array manifold. The composite array manifold (c(9)) can be produced with a 
Hadamard product (i.e., element-wise product) of the two characteristics of the array and 
expressed as: 

c(9) = a(9)©b(9) = [al (6)1)1(9) a2(0)b2(9) . . . aN(9)bN(9)]T (3) 
Where © is the Hadamard product, a(9) is the vector-valued collection of the phase 
differences of equation (1) and b(9) is the vector- valued collection of the antenna element 
directivities for an antenna array of equation (2). While the foregoing composite array 
manifold production is used in this detailed description of the drawings, any number of 
techniques can be used in accordance with the present invention to exploit the different 
features of array manifolds for DOA estimation. 

With particular reference to the enlarged view 52 of FIG. 3, the differential electrical 
phase (4>) of the signal from the remote unit is a function of the DOA (9), separation 
distance (d) 54 between the first antenna 34 and second antenna 36 forming one of the 
antenna pair (i.e., inter-element spacing), and wavelength (X) of the signal as follows: 

<|>(9) = <)>2(9)-<|)l (9) = 27i(d/X)sin(9) (4) 
While a separation distance (d) 54 that is approximately less than one-half the wavelength 
(k) of the signal frequency produces a unique mapping between the differential electrical 
phase (<j)) and the DOA (9) within a one hundred and eighty degree (180°) range, a 
separation distance (d) 54 that is approximately greater than or equal to one-half of the 
wavelength (X) (i.e., irregular antenna element spacing) introduces ambiguities in the 
mapping of the differential electrical phase (<|0 observed by the receiving antenna elements 
to the DOA (9). 

For example, the relationship between the differential electrical phase (<|0 56 and the 
DOA (9) 58 is shown in FIG. 5 for a two-element array with two-wavelength inter-element 
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spacing (d) and phase wrapping for the differential electrical phase (<j>) 56 having a 
magnitude greater than one hundred and eighty degrees (180°). For this two-wavelength 
inter-element spacing (d), there are four possible DOA solutions (60,62,64,66) for each of 
the differential electrical phases (((>) 56. Therefore, to identify a DOA estimate, the 
apparatus for determining a DOA using these phase differences preferably distinguishes 
between the multiple solutions (e.g., distinguishes between the four DOA solutions 
(60,62,64,66)). Most preferably, the apparatus for determining the DOA of the RF signal 
distinguishes between multiple solutions that have the least amount of DOA separation. 
For example, fifteen degrees (15°) of ambiguity is preferably resolved for the two 
wavelength inter-element spacing, as the least amount of separation between two of the 
multiple solutions is thirty degrees (30°). (See FIG. 6 for examples of the resulting least 
DOA separation between multiple solutions for various inter-element spacing). 

Referring to FIG. 7, the sectorized base site 24 of FIG. 1 is shown in additional 
detail that is configured to identify a DOA estimate of a remote unit with an antenna array 
having irregularly spaced antenna elements. While the following detailed description of the 
drawings is directed to a transceiver having the DOA apparatus, it should be understood that 
the DOA apparatus can be a separate component of the sectorized base site or located 
outside of the sectorized base site, including with the remote unit. Furthermore, while the 
following detailed description of the drawings is directed to identifying the DOA with the 
first antenna and second antenna in order to maintain simplicity and clarity, it should be 
understood that the description is applicable to other antennas of the antenna array and any 
combination thereof. In addition, while the following detailed description of the drawings is 
described with reference to a phased array antenna having irregular element spacing, the 
present invention is applicable to antenna arrays having regular element spacing (i.e., there 
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is a separation distance (d) that is approximately less than one-half of the wavelength (X) of 
the RF signal). 

The RF signal 32 transmitted by the remote unit 22 is received at the first antenna 34 
and the second antenna 36 within the sectorized base site 24. The first antenna 34 and 

5 second antenna 36 are configured to intercept electromagnetic radiation of the RF signal 32 
within a selected frequency band. The first antenna 34 and second antenna 36 can be any 
number of electromagnetic interception devices, including, but not limited to a wave guide, 
a coaxial cable, an optical fiber or an infrared frequency transducer. The electromagnetic 
radiation of the RF signal 32 that is intercepted by the first antenna 34 and the second 

1 0 antenna 36 is converted to a first analog signal 68 with a first RF converter 70 and a second 
analog signal 72 with a second RF converter 74. The first analog signal 68 and second 
analog signal 72 are digitized with a first analog-to-digital (A/D) converter 76 and second 
A/D converter 78, respectively, to produce a first digitized signal 80 and a second digitized 
signal 82 representing the electromagnetic radiation intercepted by the first antenna 34 and 

15 the second antenna 36 within the intercepted frequency band. Alternatively, each of the 
antennas (34,36) may be coupled to a single RF converter and A/D converter for conversion 
and digitization. 

The RF signal 32 received by the first antenna 34 and the second antenna 36, 
converted by the first RF converter 70 and second RF converters 74, and digitized with the 

20 first A/D converter 76 and second A/D converter 78 includes multiple scattered rays 
produced by multi-path scattering of the RF signal 32 during transmission from the remote 
unit, with each ray of the multiple scattered rays being a separate multi-path reflection of the 
RF signal 32. The multiple scattered rays produced by the multi-path scattering are 
generally arriving from different directions, at different time delays, and generally have 

25 different complex valued amplitudes (i.e., gain and phase components). Therefore the RF 
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signal 32 received from the remote unit 22 generally includes multiple scattered rays having 
different complex amplitudes, times of arrival, and directions of arrival. The collection of 
voltages produced by each scattered ray (ri) on the antenna array elements can be expressed 
as: 

ri(t) = c(9i)s(t-Ti) (5) 
Where xi is the time delay of ray i, 9i is the DOA of ray i, s(t) is the transmitted signal and 
c(0i) is the composite array manifold of equation (3). A model for the vector-valued 
received signal r(t) with multi-path scattering incorporating M rays can be expressed as: 

r(t) = rl (t) + r2(t) + . . . + rM(t) (6) 
In a preferred exemplary embodiment of the present invention, the RF signal 32 is an 
encoded spread-spectrum digital signal having a multiplicity of frequency and time 
overlapping coded signals from the remote unit and other remote units within the 
communication system. The multiplicity of frequency and time overlapping coded signals 
are typically transmitted simultaneously at the same or substantially similar radio frequency 
and are distinguishable through specific modulation and spreading. In other words, the RF 
signal 32 received at the sectorized base site 24 is a composite signal of individual signals 
transmitted within the communication system, and the signal of an individual remote unit is 
generally distinguishable only after despreading and demodulation. Therefore, the first 
digitized signal 80 and second digitized signal 82 are presented to a first despreader 84 and 
first demodulator 86 and a second despreader 88 and second demodulator 90, respectively. 

More specifically, the first and second digitized signals (80,82) are provided to the 
first and second despreaders (84,88) for despreading and generation of a first despread 
signal 92 and a second despread signal 94. The first despread signal 92 and second 
despread signal 94, representing the transmission from an individual remote unit, are then 
demodulated with the first demodulator 86 and second demodulator 90 into a first 

10 
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demodulated signal 96 and a second demodulated signal 98 having in-phase (I) and 
quadrature phase (Q) components. 

As previously discussed in this detailed description of the drawings, the RF signal 32 
received from the remote unit 22 generally includes multiple scattered rays having different 
amplitudes, directions of arrival and time delays, therefore the in-phase (I) and quadrature 
phase (Q) components of the first and second demodulated signals (96,98) have also 
undergone multiple reflections that result in echoes of the in-phase (I) and quadrature phase 
(Q) components. Referring to FIG. 8, examples of the in-phase (I) and quadrature (Q) 
components of the first and second demodulated signals (96,98) are shown in greater detail. 
The echoes of the in-phase (I) and quadrature phase (Q) components are generally of 
different amplitudes and different time delays, resulting in multiple scattered rays for each 
in-phase (I) and quadrature phase (Q) component of an individual remote unit. 

More particularly, a first prompt ray 100 and a first plurality of echo rays 102 result 
for each in-phase (I) and quadrature phase (Q) component of the first demodulated signal 96 
and a second prompt ray 104 and second plurality of echo rays 106 result for each in-phase 
(I) and quadrature phase (Q) component of the second demodulated signal 98. The first and 
second prompt rays (100,104) and first and second plurality of echo rays (102,106) 
generally have different amplitudes, DOAs and time delays. In a preferred embodiment, the 
despreading and demodulation functions isolate the plurality of time-delayed rays received 
at the antenna array. A matrix-valued signal (R) results from the despreading and 
demodulation, contains representations of the scattered plurality of rays (100,102,104,106) 
of the first and second demodulated signals (96,98) and can be expressed as: 

R = [ c(01)h(Tl) c(92)h(x2) c(03)h(x3) ] = [rl r2 r3] (7) 
The first column (rl) of the matrix-valued signal (R) represents the plurality of prompt rays 
(100,104) received by the antenna array. The subsequent columns (r2,r3,r4...) represent the 

11 
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subsequent time-delayed echoed rays (102, 106) as received by the antenna array. The 
scaling values (h(xi)) represent a scaling due to the time delay of the echoes that is common 
across substantially all antenna elements. The values contained in the matrix-valued signal 
(R) are provided to a signal combiner 108, decoder 110 and DOA computer 112 for 
5 subsequent processing. 

In a preferred exemplary embodiment of the present invention, the signal combiner 
108 and decoder 1 10 are configured to perform CDMA combining and decoding of in-phase 
(I) and quadrature phase (Q) components of each remote unit as described in the IS-95A 

Q Reference Document. For example, the combiner 108 receives the in-phase (I) and 

P 

Hi 10 quadrature phase (Q) components of the first and second prompt rays and the first and 

f'\ second plurality of echo rays and combines these scattered plurality of rays to form one 

T coherent in-phase (I) and quadrature (Q) signal for each remote unit. This coherent in-phase 

p (I) and quadrature (Q) signal generated from the scattered plurality of rays 

fy 

=G (100,102,104,106) is grouped into predetermined length groups (e.g., sixty-four sample 

o 

^ 1 5 length groups) of sampled signals that are independently provided to the decoder 1 1 0, which 
is most preferably an orthogonal decoder, for subsequent decoding. In addition to 
generating the coherent in-phase (I) and quadrature (Q) signal from the scattered plurality of 
rays (100,102,104,106) for decoding with the decoder 110, the DOA computer 112 is 
configured to generate a DOA estimate from the scattered plurality of rays 
20 (100,102,104,106). 

Referring to FIG. 9, the DOA computer 112 is shown in greater detail, which is 
configured to estimate a DOA. The DOA computer 112 includes a ray selector 114, 
difference calculator 116 and angle estimator 118. The ray selector 114 receives the 
demodulated signals (96,98) having the matrix-valued signals (R) from the demodulators 
25 and selects the prompt rays (rl=[rl 1 rl2]T) of the matrix- valued signal (R) or one of the 

12 
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subsequent time-delayed echoed rays (ri=[ril ri2]T) for DOA estimation. The selected rays 
(ri) 100 are provided to the difference calculator 116 for determination of an amplitude 
difference 120 and a phase difference 122 between the values of the selected rays (ri) 100 
for the two antenna elements (e.g., ril and ri2). 
5 Referring to FIG. 10, the difference calculator 116 is shown in greater detail. The 

first value (ril) of the selected rays 100 is provided to a first magnitude calculator 124 and a 
first phase calculator 126 and the second value (ri2) of the selected rays 100 is provided to a 
second magnitude calculator 128 and a second phase calculator 130. The first magnitude 
calculator 124 and second magnitude calculator 128 compute the magnitude of the first 
10 value (|ril|2) 132 and the magnitude of the second value (|ri2|2) 134, respectively. The 
magnitude of the first value (|ril|2) 132 and magnitude of the second value (|ri2|2) 134 are 
provided to a magnitude difference calculator 136 that produces the amplitude difference 
120 between the first value (ril) and the second value (ri2) (i.e., (|ril|2/|ri2|2) or (10 
logl 0(|ril |2)-1 Ologl 0(|ri2|2))). In addition, the first phase calculator 126 and second phase 
1 5 calculator 1 30 compute the phase of the first value (Zri 1 ) 1 40 and phase of the second value 
(Zri2) 142, respectively, which are provided to a phase difference calculator 144 that 
produces the phase difference 122 between the first value (ril) and the second value (ri2). 
The phase difference 122 is provided to the angle estimator 118 for calculation of multiple 
DOA estimates and the amplitude difference 120 is also provided to the angle estimator 118 
20 for selecting one of the multiple estimates as shown in FIG. 9. 

Referring to FIG. 11, the angle estimator 118 is shown in greater detail. The 
magnitude difference is passed to a first DOA solution estimator 124, which is used to 
compute a first DOA estimate 126 (61'). The first DOA estimate 126 (01') can be 
computed with any number of techniques, including the methods described in U.S. Patent 
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5,786,791, issued to Eugene J. Bruckert on July 28, 1998 and assigned to Motorola, Inc., 
which is hereby incorporated by reference. 

The first DOA estimate 126 (91') and antenna array element directivity data 128, 
corresponding to the antenna array element directivities (b(9)) discussed with reference to 
equation (2), are provided to a phase calibration computer 130. The phase calibration 
computer 130 computes a phase calibration 132, which is used to remove the effect of the 
antenna element directivities 128 on the phase difference 122, delivering a calibrated phase 
difference 136. 

The phase calibration 132 can be calculated using any number of techniques. For 
example, as can be appreciated by one of ordinary skill in the art, the selected column (ri) 
100 for the two antenna elements can be expressed in greater detail as follows: 

ri = [cl(0i)c2(ei)]Th(xi) (8) 
or ri = [al (Gi)b 1 (Gi) a2(0i)b2(9i)] T h(xi) (9) 
Where bl(9i) and b2(9i) are the antenna element directivity of the first antenna element and 
second antenna element, respectively (i.e., the antenna array element directivity data 128), 
h(xi) is a complex amplitude for the selected time delayed echo which is common to the 
antenna elements, and ak(9i) represents the complex valued phase shift observed on antenna 
k due to the positioning of antenna element k and the DOA 9i, as previously discussed in 
this detailed description of the drawings. In a preferred embodiment, the phase difference 
calibration 132 ((j)') is the phase difference between the directivities of the antenna array 
elements measured at the first DOA estimate 126 (01 '), which can be expressed as follows: 

4)' = Z(b2(91')/bl(91')) (10) 
Where b2(91')and bl(91') are the antenna element directivities of the first antenna element 
and second antenna element at the first DOA estimate 126 (91'), respectively. 

14 
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The phase difference calibration 132 (<]>') is provided to the phase difference 
calibrator 134, which preferably computes the calibrated phase difference 136 (cp") from the 
phase difference 122 (cp') and the phase difference calibration 132 (<[>') as follows: 

(p"=(<p' - <j)' ) mod 2% (11) 
The calibrated phase difference 136 (cp") is provided to an ambiguous DOA solution 
generator 138, which is also configured to receive antenna array element geometry data 138 
that is sufficiently related to the antenna array geometry for use in deriving the vector- 
valued collection of the phase differences (a(9)). 

The ambiguous DOA solution generator 140 calculates K ambiguous DOA values 
142 (e"=[61" 92" 9K"]T), preferably calculates the majority of ambiguous DOA values, 
and most preferably calculates substantially all or all the ambiguous DOA values, which 
yield the observed calibrated phase difference 136 (cp") for the given antenna array element 
geometry data 138. For example and with reference to FIG. 5, the multiple ambiguous 
DOA values 142 (i.e., the K ambiguous DOA values) can be selected with a reverse table- 
lookup operation of the ambiguous DOA values 142 (9"=[91" 92" 9K"]T) along the DOA 
axis 59 that yields the observed calibrated differential electrical phase along calibrated 
differential electrical phase axis 57. Referring to FIG. 11, the ambiguous DOA values 
(9"=[91" 92" 9K"]T) 142 and the first DOA estimate 126 (91') are provided to the final 
DOA estimator 144. 

The final DOA estimator 144 selects a final DOA estimate 146 from the ambiguous 
DOA values 142 using the first DOA estimate 126. The final DOA estimate 146 is 
preferably selected from the ambiguous DOA values 142 that have a value with the least 
numerical difference from the first DOA estimate 126. Therefore, the final DOA estimator 
144 utilizes the antenna array element directivities and the antenna array geometry to deliver 

15 
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the final DOA estimate 146. 

While the invention has been particularly shown and described with reference to a 
particular preferred exemplary embodiment of the present invention, it will be understood 
by those skilled in the art that various changes in form and details may be made therein 
5 without departing from the spirit and scope of the invention, and it is intended that all such 
changes come within the scope of the following claims. 
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